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Experimental measurements of angular correlations in three-jet final states have recently been 
performed in DIS. Next-to-leading-order QCD results for these observables are affected by large 
theoretical uncertainties in the kinematic region of the data. We discuss the effects of multiple 
QCD radiation at higher order using parton-shower methods based on transverse-momentum de- 
pendent parton distributions and matrix elements. We present Monte-Carlo results for azimuthal 
two-jet and three-jet distributions, and discuss the comparison with experimental data. 
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Three-jet final states 



The study of hadronic final states with multiple jets at the Tevatron, HERA and LHC colliders 
relies both on perturbative multi-jet calculations (see ^ for a recent overview) and on realistic 
event simulation by parton-shower Monte-Carlo generators (see e.g. [^, ^]). These complex pro- 
cesses, characterized by multiple hard scales, are potentially sensitive to effects of QCD initial-state 
radiation that depend on the finite transverse-momentum tail of partonic matrix elements and dis- 
tributions [^, 0]. In perturbative multi-jet calculations truncated to fixed order in a, \^, finite-k^ 
contributions are taken into account partially, order-by-order, through higher-loop corrections. On 
the other hand, in standard shower Monte-Carlos such as Herwig [||] and Pythia [^], based on 
collinear evolution of the initial-state jet, finite-kj^ contributions are not included, but rather cor- 
respond to corrections [^, ^ to the angular or transverse-momentum ordering implemented in the 
parton-branching algorithms. 

In this article we discuss the role of initial-state radiative effects in the study of angular cor- 
relations for multi-jet processes, concentrating on the case of three-jet DIS production, for which 
new experimental data have recently appeared [^] and next-to-leading-order (NLO) calculations are 
available [ 10 ] . In the first part we recall the experimental results [^ and observe that while inclusive 
jet cross sections are reliably predicted by NLO perturbation theory, jet correlations are affected 
by significant higher-order corrections, increasing as the momentum fraction x decreases, giving 
large theoretical uncertainties at NLO. Sizeable contributions arise from regions [^ with three 
well-separated hard jets in which the partonic lines along the decay chain in the initial state are not 
ordered in transverse momentum. In the second part we move on to the description of multi-jet 
correlations by parton showers that incorporate finite-k^ corrections to the transverse-momentum 
ordering [|n|]. The results compare well with experimental data, and give quite distinctive features 
of the associated jet distributions compared to standard showers such as HERWIG. 



1. Jet distributions from NLO calculations 



DIS multi-jet distributions associated with Q > 10 GeV and 10^ <x < 10^" have recently 
been measured by the ZEUS collaboration [Q], and compared with next-to-leading-order calcula- 
tions [10|. Results for di-jet distributions are shown in Fig. |l] [^. The plot on the left in Fig. |] 
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Figure 1: (left) Bjorken-x dependence and (right) azimuth dependence of di-jet distributions at HERA as 
measured by ZEUS [||]. 



shows the x-dependence of the di-jet distribution integrated over A0 < 27r/3, where A0 is the az- 
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imuthal separation between the two high-^T- jets. The plot on the right shows the di-jet distribution 
in A(^ for different bins of x. The overall agreement of data with NLO results is within errors. 
However, the variation of the predictions from order-a^^ to order-a^^ is significant. Given the large 
difference between order-a^^ and order-a^^ results, the theoretical uncertainty at NLO appears to be 
underestimated by the error band obtained from the conventional method of varying the renormal- 
ization/factorization scale. ^ 

The stability of predictions for the jet observables in Fig. |l| depends on a number of physical 
processes. Part of these concern nonperturbative dynamics, including jet clustering and hadroniza- 
tion^, part concern radiative corrections at higher order. Fixed-order calculations beyond NLO are 
not within present reach for multi-jet processes in hadronic collisions. Resummed calculations of 
contributions from multiple infrared emissions are performed in [p^]. These contributions are en- 
hanced in the region where the two high-^r jets are nearly back-to-back, and are also taken into 
account by parton-shower methods as in Herwig [||]. Note however that sizeable corrections in 
Fig. [I] arise for decreasing A0, where the two jets are not close to back-to-back and one has ef- 
fectively three well-separated hard jets |H]. The corrections increase as x decreases. By analyzing 



the angular distribution of the third jet, Ref . [11] finds significant contributions from configura 



tions where the transverse momenta in the initial-state shower are not ordered. These contributions 
are not included in standard parton showers, e.g. HERWIG |^ and Pythia [^]. The next section 
describes results of parton-shower calculations that take these contributions into account. 

2. Parton showers with k^ -branching and jet correlations 

Corrections to the initial-state shower due to non-ordering in k^ can be incorporated in Monte- 
Carlo event generators [|^, 14, |l^, [T^, with logarithmic accuracy for x <C 1, by implementing 



unintegrated parton distributions defined through high-energy factorization []8]]. Although none of 
the above generators are as developed as standard Monte-Carlos like HERWIG [Q] or PYTHIA |^], 
they have the potential advantage of a more accurate treatment of the initial-state radiative effects 
that the observations of Sec. [I] suggest to be relevant for multi-jet correlations. 

Fig. H shows results [jlT]] for the azimuthal distribution of di-jet and three-jet cross sections 
obtained by the kj^-shower Monte-Carlo CASCADE [|l^] and by HERWIG ^, compared with the 
measurement |^. The shape of the distribution is different for the two Monte-Carlos. CASCADE 
gives large differences from HERWIG in the region where the azimuthal separations A<p between 
the leading jets are small. It becomes closer to HERWIG as A<p increases. This is consistent with the 
expectation that both Monte-Carlos give reasonable approximations near the back-to-back region. 
The description of the measurement by CASCADE is good, whereas HERWIG is not sufficient to 
describe the measurement in the small A(p region. 

In the kj^-shower calculation both the parton distributions (unintegrated pdf 's, fitted from ex- 
periment) and the hard matrix elements (ME's, computed perturbatively) are transverse-momentum 
dependent. The parton branching includes regions that are not ordered in kj^ along the initial-state 
decay chain. Fig. ^ shows different approximations to the azimuthal dijet distribution normalized 



'Besides angular correlations, sizeable uncertainties at NLO also affect other associated distributions such as mo- 
mentum correlations. On the other hand, NLO results are much more stable for inclusive jet cross sections [H]. 
See ni in and references therein for more discussion of these effects. 
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Figure 2: Angular jet coiTelations [yj obtained by CASCADE and Herwig, compai-ed with ZEUS data | 
(left) di-jet cross section; (right) three-jet cross section. 



to the back-to-back cross section. The sohd red curve is the full result. The dashed blue curve is 
obtained from the same unintegrated pdf 's but not including the transverse-momentum dependence 
of the hard ME. We see that the high-k^ component in the hard ME is essential to describe jet cor- 
relations particularly for small A^. The dotted violet curve is obtained from the unintegrated pdf 
without any resolved branching, corresponding to nonperturbative, predominantly low-kj^ modes. 
The results of Fig. ^illustrate that the kj^ -dependence in the unintegrated pdf alone is not sufficient 
to describe jet production quantitatively. 
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Figure 3: Azimuthal distribution normalized to the back-to-back cross section [|ri|]: (solid red) full result 
(u-pdf © ME); (dashed blue) no finite-kx correction in ME (u-pdf © MEcjiiw.)'^ (dotted violet) u-pdf with 
no resolved branching. 



Jet multiplicities and jet momentum correlations are analyzed in [ ]11[ ] along lines similar to 
those discussed above for azimuthal distributions. The largest differences between HERWIG and the 
kx -showering are found in the contribution of high multiplicities and correlations in the transverse- 
momentum imbalance between the leading jets. 

In summary, the results give a consistent paiton-shower picture of distributions associated 
to multi-jet production, including correlations. This is expressed in terms of unintegrated pdf 's 
convoluted with transverse-momentum dependent hard kernels.^ The physical picture reflects the 



The u-pdf s are defined gauge-invariantly [g] for small x. General operator definitions, including the region x ~ 1, 
and issues of lightcone singularities are discussed e.g. in jfo, nM and references therein. 
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growth of the kj^ transmitted along the initial-state jet, and includes corrections to the collinear 
ordering implemented in standard showers such as Herwig and Pythia. We have considered 
production of jets in the non-forward region. Results for jet correlations are less dependent on 
details of u-pdf evolution models than in the case of forward-region observables. Owing to the 
large phase space available for jet production, the DIS kinematic region under consideration may 
be relevant for understanding the extrapolation of initial-state showering effects to the LHC, despite 
the lower energy. 
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